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D Abstract 

We analyse the neutral scalar sector of the MSSM without R-parity. Our analysis is performed 
for a one-generation model in terms of "basis- independent" parameters, and includes one-loop 
^ \ corrections due to large yukawa couplings. We concentrate on the consequences of large Rp 

d ' violating masses in the soft sector, which mix the Higgses with the sleptons, because these are 

only constrained by their one-loop contributions to neutrino masses. We focus on the effect of 
i?p-violation on the Higgs mass and branching ratios. We find that the experimental lower bound 
on the lightest CP-even Higgs in this model can be lower than in the MSSM. 
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1 Introduction 



Supersymmetry(SUSY) |T|, ^ |] is a popular extension of the Standard Model (SM), that introduces 
new scalar partners for SM fermions and new fermionic partners for SM bosons. A consequence of 
the enlarged particle content of SUSY models is that baryon (B) and lepton (L) number are not 
automatically conserved in the renormalisable Lagrangian. In the Standard Model, gauge invariance 
implies that B and L are conserved in any terms of dimension < 4; this is no longer the case in SUSY, 
so a discrete symmetry is often imposed to forbid the unwanted interactions that violate B and/or L. 

There are a variety of discrete symmetries that can be imposed to remove the renormalisable 
B and L violating terms from the SUSY Lagrangian. The most common is i?-parity 0, under which 
particles have the charge Rp = (_ i)3-B+i+25^ where 5* is the spin. SM particles are even under this 
transformation, and SUSY partners are odd, which forces SUSY particles to always be made in pairs 
and forbids the Lightest Supersymmetric Particle (LSP) from decaying. 

Alternatively, one can allow the B and L violating interactions to remain in the SUSY Lagrangian, 
and constrain the couplings to be consistent with present experimental data. The renormalisable Rp 
violating couplings violate either B, or L. If both types of coupling are simultaneously present, they 
can mediate proton decay, and are therefore constrained to be very small 0. So in this paper, we 
will assume that the B violating couplings are absent — forbidden by some other symmetry — and only 
consider the L violating couphngs. These are particularly interesting, because L violation is observed 
in neutrino masses. 

The renormalisable Rp violating interactions have a variety of phenomenological consequences . 
These include generating majorana neutrino masses, mediating various flavour and lepton number 



violating processes |T0|, |Tl|], and modifying the signatures of supersymmetric particles at colliders 
11, ^ In particular it allows the lightest supersymmetric particle (LSP) to decay 0, |l5|. It can also 
modify the Higgs sector. 

The Higgs sector of the i?p-conserving MSSM has been extensively studied |16|, [l^, |l^, |19|, 
2T| , [2^ , with a lot of emphasis on both one-loop [^, ^ |2^, ^ and more recently on two- loop 



effects [19|, ^7\, ^ |2^, ^ ^ to the lightest Higgs boson mass. The most relevant one-loop effects 
due to the large top-quark Yukawa coupling are from the stop-top sector. There are several different 
approaches that have been utilised to incorporate these loop effects: effective potential methods, 
renormalisation group running, explicit diagrammatic calculations (see e.g. for a review). The 
effective potential, which we use here, can in a simple way take into account the most relevant effects 
although it does not incorporate any momentum-dependent contributions. 

A Higgs boson could be the next particle discovered at accelerators. It is therefore interesting 
to study its properties in various extensions of the Standard Model, in particular SUSY. One of the 
advantages of the supersymmetric Standard Model for cosmology is that baryogenesis may be possible 
at the electroweak phase transition — if the Higgs is light enough []39| , ^ . However, as the experimental 
lower limit on the Higgs mass increases, the parameter space remaining in the MSSM for baryogenesis 
is reduced. Adding Rp violation can decrease the experimental lower limit on the Higgs mass, which 
could increase the available parameter space for electroweak baryogenesis. 

In this paper, we study the neutral Higgs sector at one-loop in the _Rp-violating MSSM with 
one generation of quarks and leptons, since this toy model already contains the main effects of the 
complete three generation case. We vary bilinear and trilinear Rp violating parameters over their 
experimentally allowed ranges, and discuss how this can change the masses of the neutral CP-even 
scalar bosons and the branching ratios of the lightest one, hi. The violating Higgs sector has been 
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studied by numerous authors: novel decays of both neutral |^ ^ and charged scalar bosons have 
been analysed in the context of bi-linear Rp violation, and the mass matrices of the Higgs sector have 
been derived, considering only the effect of bi-linear terms |37] and in the general case with both bl- 
and tri-linear couplings |3^. Our analysis differs from previous treatments in that we include one-loop 
yukawa corrections to the Higgs masses, and we parametrise Rp violation in a basis-independent way 
that avoids much possible confusion about what is a lepton/slepton in a lepton number non-conserving 
theory. 

The next section of this paper introduces our notation and discusses the basis-independent approach 
to Rp violation. The third section is devoted to experimental constraints on the Rp violating parameters 
in our model, largely from neutrino masses. In the fourth and fifth sections, we calculate the masses 
and various branching ratios for the CP-even Higgses at one loop. We present our results in section 
six. The first appendix contains the one-loop Higgs mass matrices in an arbitrary basis. The second 
appendix contains the same information, but in the basis where the sneutrino vev is zero (to one loop). 
The third appendix contains a few useful but long formulae. 



2 Basis dependence of the Lagrangian 

In the SM, the Higgs and leptons have the same gauge quantum numbers. However, they cannot mix 
because the Higgs is a boson and the leptons are fermions. In a super symmetric model this distinction 
is removed, so the down-type Higgs and sleptons can be assembled in a vector Lj = {Hd, Li) with 
J : O..Ng = the number of generations. We write vectors in Lj space with a capitalised index J or as 
vectors v, and we write matrices in Lj space in bold face m. Using this notation, the superpotential 
for the supersymmetric SM with Rp violation can be written as 

W = fi'^H^Lj + X'/'LjLkEI + Xi^'^LjQpDl + hf'H^QpU'^ (1) 

The Rp violating and conserving coupling constants have been assembled into vectors and matrices in 
Lj space: we call the usual /i parameter /xo, and identify the usual ei = ^i, = X^J^, X^^^ = X^J'', 

~ ■^'b"^} and A *^"^ = A^*^^. Lower case roman indices k and p, q are lepton and quark generation 
indices. In the body of the paper, we will work in a one generation model, so = A°^, hb = A°, 
and a' = A^ and now the capitalised indices run from 0..1, and 1 corresponds to the third lepton 
generation. We often write d and L (for down- type Higgs and slepton) rather than and 1. Q, U'^ 
and D'^ are the third generation quark superfields. In the one-generation model, there is no XLLE'^ 
interaction (because A is antisymmetric on the capitalised indices). 

We also include possible Rp violating couplings among the soft SUSY breaking parameters, which 
can be written as 

Vsoft = -fHlH^ + -L'^[ml]jKL'' + B'H^Lj + AtH^QU'' + AiLjQD' + Ai'^LjLKE'^ + h.c. (2) 

Note that we have absorbed the superpotential parameters into the A and B terms; e.g. we write 
B^HuHd not B^ ^^HuHd Q. We abusively use capitals for superfields (as in (|^)) and for their scalar 
components. 

^We do this because Bj is a vector — a one index object — in {Lj} space. ^From this perspective, giving it two indices 
can lead to confusion. 
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The reason we have put the Higgs Hd into a vector with the sleptons, and combined the Rp - 
violating with the Rp conserving couphngs, is that the lepton number violation can be moved around 
the Lagrangian by judiciously choosing which linear combination of hypercharge = -1 doublets to 
identify as the Higgs/higgsino, with the remaining doublets being sleptons/leptons. This can create 
some confusion when one tries to set experimental constraints on lepton number violating couplings; 
it makes little sense to set an upper bound on a coupling constant that can be made zero by a basis 
rotation. 

If one calculates a physical observable as a function of measurable quantities, then the basis in which 
one does the intermediate steps of the calculation is irrelevant. However, if one computes observables 
as a function of Lagrangian quantities, as is common in Supersymmetry, it can be important to specify 
the basis chosen in the Lagrangian. In SUSY theories with lepton number violation, there are various 
possible choices for what one identifies as a lepton/slepton in the Lagrangian, and the interactions 
that are "lepton number violating" depend on this identification. However, this freedom to redefine 
what violates L is deceptive, because phenomenologically we know that the leptons are the mass 
eigenstate e, /i and r , so we know what lepton number violation is. There are two possible approaches 
to this fictitious freedom; either one chooses to work in a Lagrangian basis that corresponds to the 
mass eigenstate basis of the leptons, or one can construct combinations of coupling constants that are 



independent of the basis choice to parametrise the Rp violation in the Lagrangian [|12], ^ 
These invariant measures of Rp violation in the Lagrangian are analogous to Jarlskog invariants which 
parametrise CP violation. 

The standard option is to work in a basis that corresponds approximately to the mass eigenstate 
basis of the leptons. For instance, if one chooses the Higgs direction in Lj space to be parallel to /ij, 
then the additional bilinears in the superpotential Hi will be zero. In this basis, the sneutrino vevs are 
constrained to be small by the neutrino masses, so this is approximately the lepton mass eigenstate 
basis. Lepton number violation among the fermion tree-level masses in this basis is small by construc- 
tion, so it makes sense to neglect the bilinear Rp violation, or treat the small Rp violating masses as 
"interactions" within perturbation theory, and set constraints on the trilinears, as is commonly done 
(for a review, see .e.g. |T0|. For a careful analysis including the bilinears, see p5|). 



In this paper, we present our results in terms of basis-independent "invariants". We also give 
explicit results in the basis where the sneutrino does not have a vev, which is close to the lepton 
mass eigenstate basis. This is to present our calculation in a familiar way. The advantage of the 
first approach is that we can express Higgs masses and branching ratios in terms of inputs that are 
independent of the choice of basis in the Lagrangian. The drawback is that the "invariants" can 
appear unwieldy and forbiddingly complicated. However, since we work in a model with only one 
lepton generation, the linear algebra is tractable. 

The aim of the "basis-independent" approach is to construct combinations of coupling constants 
that are invariant under rotations in Lj space, in terms of which one can express physical observables. 
By judiciously combining coupling constants one can find "invariants" which are zero ii Rp is conserved, 
so these invariants parametrise Rp violation in a basis-independent way. For instance, consider the 
superpotential of equation (|l|) in the one generation limit, I : 0..1. It appears to have two Rp violating 
interactions: fiiH^L and X'LQV^. It is well known that one of these can be rotated into the other by 
mixing and L [Q. If 
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L' = I --Hd , ^° L , (3) 

then the Lagrangian expressed in terms of H'^ and L' contains no HuL' term. One could instead dispose 
of the X'LQD^ term. The couphng constant combination that is invariant under basis redefinitions in 
{Hd, L) space, zero if R parity is conserved, and non-zero if it is not is /iqA' — /i^/ii = (/io, /ii) A {hd, A'). 

In this paper, we are interested in Rp violating effects in the Higgs sector, so we are interested in 
constructing invariants involving Bj, the L j mass matrix [m|] ji^ = [mWjK + jjJ'K, and the L j vev 
vj = (i^j). The vev vj is a dependent variable, fixed by Bj and [m|] jj^ in the minimisation conditions. 
In an arbitrary basis, there are therefore two R^ violating masses in the Higgs sector: Bx and [m|^]oi. 
However one can always choose the basis such that one of these parameters is zero, so we expect only 
one independent invariant parametrising i?p violation in the (tree-level) Higgs mass matrices. 

There is R^ violation in the Higgs sector if B, [m|], and v disagree on which direction in Lj space 
is the Higgs, or equivalently, if it is not possible to choose a basis where = Bj^ =\a^\dL = 0. -B is 
a vector that would like to be the Higgs — that is, if the basis in Lj space is chosen such that Hd oc B 
then Bd = \B\ and Bl = 0, so the mass matrix mixes H^ with Hd but not with L. [mj^]jK has two 
eigenvectors in Lj space, one of which would like to be the Higgs, and the other the slepton. v is also 
a candidate direction in Lj space to be the Higgs — the basis where Hd is the v direction is the basis 
where the sleptons do not have vevs. There is Rp violation if two of B, v and [m|^]ji^ do not agree on 
what is the Higgs direction. A convenient choice for the invariant parametrising this Rp violation at 
tree-level is 

R = v'\B\'-{v.Br ;5« = ^ , (4) 

where 6r is the normalised version of the parameter, varying from for no Rp violation to 1 for 
maximal Rp violation. As we will see from the minimisation conditions (equations |22| and ^3|), at tree 



level xB = —[mj] ■ v, where x is the vev of the up-type Higgs, so we can also write x^R = "^^"^ ' [ 
m|]^ ■ V ~ (v ■ [m|] • v)"^. R parametrises the Rp violation in the mass matrix relevant for the Higgs. 

is the sine of the angle [| between B and v (see figure m), which is clearly independent of the choice 
of basis in Lj space. 

There are many other invariants that parametrise Rp violation among other coupling constants. 
For instance, there is an additional invariant among the bilinears in one generation There 



are three possible directions in Lj space that could be identified as the Higgs: Bj, fxj and one of 
the eigenvectors of [m|]j/^. If these three vectors do not coincide, there should be two invariants 
parametrising the misalignment between the three vectors. One in the scalar sector, as constructed 
in equation (^, and an additional one involving fl. For instance, if /2 is misaligned with respect to v, 
mixing between neutrinos and neutralinos generates a tree-level neutrino mass ~ flAv = v ■ Xr ■ fJ^/lXrl- 
Invariants parametrising Rp violation between bilinears and trilinears can also be constructed. Since 
the upper bound on neutrino masses constrains fl A v to he small, we neglect it in this paper, and 
concentrate on the effects of 6^. 

Up to this point, we have discussed the construction of invariants using parameters from the 
Lagrangian without specifying whether they were tree-level, or computed to some loop order. We 
choose to write the invariants in terms of one-loop parameters. We do this because the invariants 
were constructed to avoid expressing measurable quantities {e.g. masses) in terms of unmeasurable 



We take the positive square root: sinry = +v7r 
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Figure 1: Non-orthogonality of the soft mass BjHuL^ and the if^-slepton vev (Lj) = vj. The angle 
T] between B and v is basis-independent. The "invariant" R = v^B'^ — {v ■ By is equal to v^B^ sin^ rj. 
The basis here is H oc v, and L cav ■ X^-. 
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basis dependent Lagrangian parameters. So we define the invariants in terms of one-loop parameters, 
because tliese are closer to wliat is pliysically measured. Tlie invariants R and 5r discussed above are 
tlierefore taken to be 

R = v'Ml-{v-M^f ■,Sr = ^, (5) 

wfiere M„ is tlie one-loop corrected version of B tliat appears in tlie CP-odd mass matrix (0). From 
tlie one-loop minimisation conditions ( p2|) and (|23|) , M„ = —'M.-v/x, where M is the one-loop version 
of m| that appears in the CP-odd mass matrix. R can therefore also be written as 

X^R = v^v-M^-v-{v-M-vf. (6) 

The drawback to using the one-loop expressions is that it is not obvious which loop corrections should 
be included. 

We will use 6r rather than R as our Rp violating parameter, because it is dimensionless and 
normalised to 1. For small this is clearly a good choice, because the magnitude of M„ is largely 
determined by its Rp conserving component (~ m a sin f3 cos/? in the MSSM). However, as 6r increases 
to 1, the magnitude of the Rp violating mass^ term |M„|a/5r can nonetheless decrease if |M„| does. 
We will see that for some parameter choices, this is the case. 

We would like to determine which are the necessary conditions on the Rp violating parameters 
to produce a substantial effect on physical observables. Hence, we do not assume in this paper that 
Bj ^ Bfij, (and [m|] ^ fh? I) as would be expected in many models of SUSY breaking. This means 
that we allow 5r to be as large as experimentally allowed. 



3 Experimental constraints 

Both low and high-energy processes can place stringent bounds (see e.g. [|, |lOl) on the -Rp- violating 



couplings which give rise to new interactions. The most relevant constraints on the Rp violating bilinear 
couplings come from neutrino masses. The trilinear A' also contributes to neutrino masses, but the 
most restrictive bound on A' comes from Z decay to hh. We now mention the contribution to neutrino 
masses due to various Rp violating parameters; the purpose of this discussion is to set bounds on our 
parameters, not to calculate the neutrino mass. 

In i?p-violating models the neutrino can acquire a mass at tree-level through mixing with the 
neutralinos and also through loops which violate lepton number by two units. In the basis where the 
sneutrino vevs are zero, the tree-level contribution can be written as ^ 



m|/ioM^ cos^ /3 .2. .^x 
'"'^^ ~ m|M^ sin 213 - M.M^f^o ^ ^ ' 

where Mi, M2 are gaugino masses. My = Mi cos^ 9w + M2 sin^ 9w, and sin^ = {flAv)/\fi\\v\. Thus the 
neutrino mass sets the constraint that ft be aligned with v, which determines the tree-level contribution, 
without imposing any constraints on the other i?p-violating parameters. 

There is also a loop contribution to the neutrino mass proportional to 6r, as discussed in |49|, |50|. 
If V and Mu are not parallel, then the Rp violation in the soft masses will mix the real (imaginary) part 
of the sneutrino with the CP-even (odd) Higgses. This introduces a mass splitting between ur = Re{i>) 
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and vi = Im{v). A neutrino mass can be generated by a neutralino-neutral scalar loop — see figure 
d^). The amplitude for this diagram is 

= ^ E '^x. (^.2 - Z^,g'/gf Y.iu ■ s.feMO, Mf, (8) 

Xj i 

We in practise neglect the sum over the four neutralinos and just include the lightest one. The Zij are 
the usual mixing angles between the neutralino mass and interaction eigenstate bases — for simplicity 
we only include the gauge coupling of the neutralino. We sum over three CP-even and two CP-odd 
scalars : Si = {hi, h2, h^, Ai, A2}. The {(z> ■ Si)} are the mixing angles between the neutrino and 
the various scalars Sj. They are basis-independent quantities which we will calculate as dot products 
in Lj space in section 5. is +1 for the three CP-even Higgses and —1 for the CP-odd. Bq is a 
Passarino-Veltman function: 

f (P'^ k 1 M ^ / \ 

'''^'^''-<^ = -''''''^^oJ i2.y^lik + qr-rnlW^ . (9) 

There are divergent and scale-dependent contributions to Bq in addition to the right hand side of 
equation (P); however these cancel in the sum over scalars Si in equation (H). 



u » ■ X ■ ■ < u 

X 



Figure 2: Contribution to the one-loop neutrino mass from bilinear Rp violation in the soft masses. 
This diagram is possible because the sneutrinos mix with the Higgses. 

The dependence of rrij, on different parameters can be understood in various limits. As 5/? ^ 
two of the CP even neutral scalars, hi and hj, become h and H of the MSSM, the third CP even 
scalar h^ becomes un, and Ai A oi the MSSM while A2 — I'l- The overlap between the neutrino 
and the MSSM Higgs {hi, hj, Ai} goes to zero (we will show in the next section that it is cx Sr), and 
z/ ■ zz/i; ~ z> ■ z// — i> 1 . The real and imaginary parts of the sneutrino contribute to the sum with opposite 



sign; we expect m^^ — mj^^ oc 6r so the sneutrino contribution will also go to zero with 6r ^ . 

The neutrino mass also decreases (for arbitrary Sfj) as either the neutralino mass or the CP-odd 
scalar mass m^a goes to infinity. We can therefore estimate 



2 



' (10) 



6471' 



4 ""X 



mv 00 



where we have assumed that as or become large, all remaining masses are of order mz- This 
is an overestimate, because it neglects cancellations in the sum (RI). For ~ 1, and most choices of 
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'^Ai , ■'^^2 5 ""^x tan/3, the neutrino mass will be < 10 MeV, so there is no bound on 6ji from the 
laboratory limit m^^ < 10 MeV. If we require rrii, < eV, as would be required by oscillation data, we 
find Sr ^ 10~^ for ~ fn^ ~ mz- 

The second type of loop diagrams involve fermion-sfermion loops. The contribution proportional 
to the trilinear coupling constant A^33 in the usual three-generation mass-eigenstate basis-analysis, can 
be expressed in a basis-invariant way as 

lOTT^ mf 

02 

where f{x) = —jz^, x = (^^^^ ' -^b is given in appendix A, and z) and A are defined in section 5. In 

the basis where the sneutrino does not have a vev, i> = (0, —1) and A = {hb, A'). Here i> is specifying 
the neutrino direction. So, A^gg will have a certain allowed upper value for a given set of the inputs 
that determine the sbottom mass parameters. 



Another bound on A^gg has been given in the literature from the calculation of Ri in which the 



allowed value of the coupling scales with right-handed soft-SUSY breaking mass m^^. For m^^ > 500 



GeV the trilinear coupling can be of order 1. Other bounds from Bq-Bq mixing or i? — > tvX [^21 

also have been studied and they also allow values of A^gg ~ 1 for sufficiently heavy right-handed 

sbottom, on the order of 300 GeV . 

The actual numerical bounds on the Rp violating coupling will depend on the input value one 
takes for the neutrino mass. If we use the experimental limit on the tau neutrino mass ~ 10 MeV, 
we can easily have a value \ ■ v ~ 1, thus its effect on the Higgs sector will be analogous to that 
of the top Yukawa coupling. In this case the bound from Ri is stronger for generic values of the im- 
parity conserving parameters. For smaller values of the neutrino mass, such that for example neutrino 
oscillation scenarios can be fulfilled, the bounds are very strong on the Rp violating couplings. 

Note that allowing A' ~ 1 in the fermion mass eigenstate basis means that A is almost perpendicular 
to V. The 5-quark mass is rrib = — (A • v)/y/2 <C 1^1, and X' = v A X/\v\. 

There are various accelerator limits on particle masses and coupling constants when i?-parity is not 



conserved (see e.g. ||TT| for a discussion). These often depend sensitively on a number of parameters, 
so are difficult to translate to the model we consider here. We will comment the LEP lower bound on 
the mass of sneutrinos with Rp violating decays in section 6. 



4 Higgs boson masses 

In this section, we calculate the Higgs boson masses using the effective potential. To do this we make 
an SU{2) rotation on the doublet, if^ — $d = eH*^ {ei2 = —1, e'^ = —1), to put the neutral 
component in the same element of the doublet as for H^. This makes it easy to compare the Rp 
conserving part of our calculation to standard two-Higgs doublet results. We also rotate the slepton 
field. So we can write 

H - ( ^^]-( ] H - ( ^'\-( + \ ii2\ 



*The bounds from references ||52|, p3|, |54| depend on whether the CKM-mixing is present in the down-quark sector. 
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(13) 



where we define x to be the vev, and Vd, vl to be the down-type Higgs and slepton vevs (in some 
arbitrary basis). We will not be concerned with the charged fields in this paper. 

From the superpotential and soft terms of equations (P and (^, the tree- level potential for the 
neutral scalar vevs is 

,2 



tree 



1^ 
+ 2" 



m, 



(14) 



where k = {g"^ + g 2)/8 ,v'^ = I^Tp = vj + vl, ml = ml + l/I^I and [mi] = [m|]jK + /ij/i^. 

We include the loop corrections due to large yukawa-type couplings, but not due to gauge couplings. 
The one-loop contribution to the potential from tops, stops, bottoms and sbottoms will be 

,2 \ 



V, 



loop 



647r2 



-12mf 





3' 






(ml\ 


3' 








+ 6mf 


In 




+ 6m| 


In 




~ 2 






~ 2 





t2 



-12mt 




+Qm' 



bi 



In 



rrir 



+ 6m 



b2 



In 



m 



b2 



(15) 



We include the bottom contributions because the Rp violating A' can be large (in the basis where the 
sneutrino does not have a vev). 

We are principally interested in the behaviour of the lightest CP-even neutral scalar — the "Higgs" . 
We would like to obtain its mass as a function of observables like the masses of the CP-odd scalars, 
and parameters like tan/5 and the "invariant" 5^ (equation ^ that parametrises Rp violation. We 
therefore need the 3x3 mass matrices for the CP-even and CP-odd Higgses at the minimum of the 
potential. 

The tree-level minimisation conditions can be written in terms of the CP-odd mass matrix elements 
([T9|) . In the absence of CP violation, the one-loop minimisation conditions expressed in terms of the 
one-loop CP-odd mass matrix have the same functional form (see equations (p2|) and (^3])). This is 
useful because it means we can impose the minimisation conditions at one loop without calculating 
either the one-loop CP-odd mass matrix or the one-loop minimisation conditions. To see this, we write 
the potential as a function of six variables: 

rrO* ttO uO* zjO t 0* t zjO zjO zjO r 



The three minimisation conditions for the potential can then be written 





dV ^^dV dCr, 







dV 



(16) 



(17) 



The CP-odd mass matrix is of the form 

j^CP-odd _ 




( Mud M^L 



Mdd MdL 
MdL Mll 



(19) 
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where the individual elements are 



' dV d^Cn 
2-^ Pin f^Aifi^i 



^1 dC„ 



(20) 



Note that our capitalised Ms have mass dimension 2. The indices i,j run from 1..3, or over u,d,L, 
and the are the imaginary parts of the scalars (see equation "J" is not an index in Lj space.) 
Second derivatives of V do not appear because they are multiplied by first derivatives of the {Cj}, 
which are zero (evaluated at 0j = 0). Since 



dx 



X- 



(21) 



(and similarly for the other derivatives of the {C„}), we see that the minimisation conditions can be 
written in terms of the CP-odd mass matrix: 



X 







Muj + 



mjKV 

X 



K 







(22) 



(23) 



We emphasize that these equations are valid in any basis, and we apply them at one-loop. 

Explicit formulae for the minimisation conditions and the mass matrix elements can be found in 
the Appendices. Appendix A contains the results for an arbitrary basis in terms of basis-invariant 
quantities. In appendix B we present the results in the basis = 0, using the familiar Lagrangian 
notation. 

The eigenvalues of the CP-odd mass matrix M are easy to obtain, since M has a zero eigenvalue. 
The two non-zero eigenvalues are 



m 



+ rr[M] ± 



(M„„ + Tr[M])2 - 4(M„„Tr[M] + det[m] - |M„| 



(24) 



In the Rp conserving limit, = nii,. When Rp is not conserved, the sneutrino as a complex field 
has Dirac and Majorana masses, so its real and imaginary parts are not degenerate. The mass of the 
imaginary part is what we identify here as niA^. By using the minimisation conditions (^2]) and (^), 
we can rewrite these masses in terms of "basis-independent" invariants (scalars in Lj space) as 



r 



+ Tr[M] ± 



^2 ^2 



Tr[M] 



4i? 



cos^ fi 



(25) 



Note that we have chosen to write and m^^ as functions of scalars in Lj space which are non- 
zero in an Rp conserving theory (such as Tr[M], v- M'-u) and scalars that are zero in an i?p- conserving 
theory (5_r). This is slightly different from choosing a basis in which one writes the masses as a part 
depending on Rp conserving couplings and a part depending on Rp violating couplings (as done for 
instance in iPHl), because for some basis choices the Rp conserving invariants depend on Rp violating 



couphngs {e.g. in the M^^ = basis, w-M-w = v\ M^rf -f Iv^vl M^^ + v\ M^l) 
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The CP-even mass matrix will be 



n=l 



where we have temporarily introduced Yi = {x,Vd,VL). Explicit formulae can be found in the Ap- 
pendices. We can express the eigenvalues of the CP-even mass matrix = {m1_^,m1^,m1,^} in terms of 
scalars in Lj space, by constructing the characteristic equation of (M' — I), and expressing the co- 
efficients in terms of invariants. We do not show the formulae (analogous to (0)), because they are too 
long to be enlightening. Another possible way to solve for {m1^} as a function of m\^,m\^,ta.ni3,6fi 
and loop corrections, is to express the matrix elements of M' in a basis-invariant way using equation 
(p6|). We plot the CP-even masses for various inputs in section 6. 

We have chosen m\^,m\^,teinP, and 6r as inputs because they are "physical". However there are 
relations between these parameters which constrain the ranges over which they can be varied. To solve 
for rrih^ as a function of our inputs, we invert equation (^) to write the basis invariant S = v ■ in 
terms of rriAi, ^R- 

S = "-^^ = + < ± ^{mX-m%y-AmXm%^) (27) 

where 7 = siia^ (36^/(1 — Sr), and the + (-) sign corresponds to > ("^Ai < "^yi2)- ^ 
cos^/3 when —>■ 0. Clearly this inner product must be a real number; to ensure that the square 
root is positive, we need 

K--iJ > SI (28) 



2mAjmyi2 sin /5 \l 1 — 6r 
so niAi and m^ij cannot be degenerate for non-zero Sr. 



5 Higgs Branching Ratios 

Including Rp violation in the Higgs sector will modify the interactions as well as the masses of the 
Higgses. Intuitively, it mixes the sneutrino with the neutral Higgses, so it can modify the amplitudes for 
Higgs production and for Rp conserving decays, as well as allowing new decay modes such as h —>■ ux^ 
and h r^^ Rp violating couplings also modify the decays of Higgs decay products. For 

instance, the LSP xo^ produced in h ^ x°z/ and h X^x" > could decay (to three fermions) within 



the detector 11^ , p!5| . It turns into a neutrino and an off-shell hi, which then decays to SM fermions. 
So if x° can be produced via an Rp violating vertex (in our case related to 6r), then it decays rapidly 
through the same vertex. 

The Higgs production and decay rates clearly cannot depend on the basis in which they are com- 
puted, so we will work in a "basis-independent" approach. We are principally interested in Rp violation 
from the scalar Higgs sector, as parametrised by the invariant of equation ^, so we will write the 
decay rates in terms of this and other invariants. There are three mass eigenstate bases in (if„, Lj) 
space that are relevant for calculating branching ratios: the CP-even mass eigenstate basis, the CP- 
odd basis, and the fermion mass eigenstate basis. Rotation angles between these bases will appear 
in the Higgs interaction vertices. We will provide expressions for these ("physical") angles which are 
independent of the basis choice in the Lagrangian. 
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In the i?p-conserving MSSM, the hghtest CP-even Higgs /i is a hnear combination of the up and 



down type neutral Higgses: h 
and an h is 



cosac 



sin a(f)f. The ZZh vertex via which LEP can produce a Z 



-(x cos a — f sin a) 



igmz 



[sm P cos a — cos /3 sin a) 



(29) 



2 cos^ 6w ' cos Q-w 

Single sneutrinos cannot be produced in the MSSM, but the Z can decay to a pair of them if kine- 
matically possible. The Z can similarly decay into a CP-even and odd Higgs, for which the vertex is 
proportional to gcos{P — a) / {2 cos 9w)- 

Adding Rp violation involving one lepton generation means the sneutrino mixes with the Higgses, 
so the lightest Higgs hi will be a linear combination of three fields hi = (cosa)(/)^ — (sin a cos 
— (sin a sin If we define the angle if with respect to the basis in Lj space where the sneutrino 

does not have a vev, then sinacos<y9 = —hi ■ v/\v\ and sin a sin (y9 = —hi ■ 0. The vector C> is the 
lepton direction orthogonal to the vev: = ■ v/\v\. U v A fl = 0, this vector corresponds to the 
charged lepton mass eigenstate |4^, which is the neutrino flavour eigenstate. We therefore call this 
direction 0. The ZZhi vertex is then a simple generalisation of (p9[) : 



— (sm p cos a + cos p- ■ hi 

cos 6w V 



and the ZhiAi vertex becomes 
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2 cos 9w 



(cos /5 cos a - {v ■ hi){v ■ Ai) - (z> • hi){0 ■ Ai)){pA - PhY 



(30) 



(31) 



where ph and pa are the momenta of the outgoing scalars. 

To evaluate the angles between the CP-even mass eigenstate basis and the zero-sneutrino-vev basis, 
we must identify the direction in Lj space corresponding to hi. The lightest eigenvector of the CP-even 
Higgs mass matrix satisfies 



ML 
ML 
ML 



ud ML 



dL 



^L ^'ll 



Ui ^ 




( Ui 


hid 




hid 


hiL ] 




\ hiL 



(32) 



where the mass matrix has primes to denote that it is the CP-even mass matrix and not the CP-odd 
matrix of equation |T^. We would like to solve this for hi = {hid, hn) Q. We can write this as two 
equations for scalars, vectors, and matrices in Lj space: 



and 



Rearranging (^4]), we find 



Muu^i + M^-hi = mlv.i 

UiM'^ + yt -hi =mlhi 
hi = ui[mlj 



(33) 

(34) 
(35) 



^ Normalised vectors wear hats, so for instance = 1. The mass eigenvectors hi,Aj are in the 3-d {Hu, Lj) space; 
hi is the projection on Lj space and — sin^ a 
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In a one generation model, this is simple to solve because the inverse of a symmetric 2x2 matrix 
N' = [mi l - M'] is N'-i = -e N'£/det(N'), where en = £22 = 0,£i2 = -S21 = -1- So 

and 

V V aet[i\ \ 

Ui = cos a can be determined from the normalisation of hi: u\ + h\ = 1. The vector = 
Mu — ra\ P sin /3 v/v+ loop corrections, and M^j = M/j + m| cos^ P^^^ + loop corrections, where 
Mu and M are from the CP-odd mass matrix ([T9|). These loop corrections, which are not presented in 
our analytic formulae, are listed in the Appendix. The loop contribution to the CP-odd mass matrix 
is implicitly included; the contribution missing from our analytic formulae is the one-loop difference 
between the CP-even and CP-odd mass matrices. Using the minimisation conditions (^21) and (^), 
we find 



u ■ hi = - — ^—prS tanpiml + m|(sin^ /9 — cos^ — — h loop corrections (3^ 
det[N \ ^ 1 — Or 



and 



V ■ hi 



-7- [{rn\_^m\^ sin P cos P — {S tan /? + m| cos /? sin /?) 



det[N' 

+m| sin P cos P{m\_^ + — S/ cos^ P)] + loops. (39) 



We do not present formulae for the loop corrections, but they are included in our numerical plots. S 
is defined in equation (p7|) ; the normalisation factor ui/det\N ] is in Appendix C. 

In the limit 5ji 0, the lightest CP-even Higgs hi can become either the MSSM Higgs h or the real 
component of the sneutrino ur. Suppose first that hi ^ h as 6r 0. Then as expected u ■ hi (x 5r. 
If hi — > vr in the (5r — > limit, then v ■ hi ^ because nih^ ^A2- v ■ hi ^ 1 m. the same limit, 
although this is less obvious because ui/det\N ] is singular. 

To calculate the contribution of the various Higgses to the neutrino mass, as discussed in the 
experimental bounds section, we need the angle mixing the neutrino with each of the Higgses: v ■ Si 
{si = {hi, h2, h^, Ai, A2}). These can be computed in the same way as z> ■ hi. For h2 and h^, the 
formulae are the same, substituting rrih^ or m/jg for m/j^. For Ai and A2, M' is replaced by M and 
by Mu in the analogue of equation (^5|). This gives 



Ai = niSm\, tan/?W — — h loops (40) 
» >' 1 — - 



R 



where the normalisation factor Ui is in Appendix C. 

There is a technical catch to this way of calculating the {0 ■ Si\ in the 5r Q limit. \i 5r = one 
of the hi, say /13, and A2 are the sneutrino so ■ A2 = v ■ h-^ = 1. This is the 5r ^ Q limit of equations 
and ( ^61) because the denominator 0, but at 5^ = the equations are singular. This can be 
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avoided by taking P ■ A2 = \J \ — {0 ■ AiY and ■ = y 1 — {0 ■ hiY — {0 ■ which follow from the 
unitarity of the rotation matrix. 

The tree-level rate for a scalar h to decay to two fermions /i and /2 through a vertex of the form 



hfl{XLPL + XRPR)f2 



where = (1 — 75) /2, is 



El - ml 



[ml -ml- ml){\\ + A^) - 4ALARmim2 



(41) 



(42) 



where E2 = (m| + m^ — mf) / {2mh) . 

The Rp violating decay rates h '^X^^'^X^ ^'^'^ both detectable if kinematically allowed, because 
X° can decay to u and an off-shell Higgs, which can then decay to SM fermions. Here we mention 
again that the neutralino/chargino is produced and decays via the same vertex which is proportional 
to 6ji. If 6r ^ but fl A V = f\, the decays h z/x°, tx~^ proceed because the mass eigenstate h 
contains a "(s)neutrino component" = D ■ h. The coupling constant for the vertex hx^u is therefore 



Ar = A 



R 



Zug'/g) ■ h 



(43) 



we can compute 

Of 



where Z diagonalises the neutralino mass matrix: ZmZ'^ = diag. Substituting in (| 
the decay rates r(/i — > x'^^) ^^id r{h —>■ x'^'^)- Note that by "/i —>■ x'^^" ^6 mean h — >■ xf^u and h —>■ x""^ 
The hibb coupling A ■ hi can be much larger in Rp non-conserving theories than in the MSSM. 
Decomposing A = (A ■ v)v/v'^ + (A ■ z>)z>, (in the (u) = basis this is A = {hb, A')), it follows that 



Ui 



detN' \ y/2mwcosl3 



[m\jn\^ sin [3 cos (3 — m\^ [S tan (3 + m| cos (3 sin (3)) 



-\-m\ sin (3 cos [3{m\^ + m\,^ — S/ cos^ (3) 
+X ■ V tan/9(m^^ + m|(sin^ /? — cos^ /5)) 



-^r\ 



+ loops 



(44) 



where [X-v] can be ~ 1 as discussed in section 3. This expression simplifies when Rp violation is small: 
\ih\ ^ h when — > and A ■ z> — > 0, then X-h\ ^ gm^j {yplmw cos /5) as expected in the MSSM. If 
/13 and A2 are the sneutrino components in the same limit, then one can check that A ■ /is 0. 



6 Results 

We express the masses and coupling constants of the CP-even Higgses in terms of tree-level input 
parameters m^i, mj^^^ tan/3 and When loop corrections are included there is an additional depen- 
dence on the soft parameters A, /i/, mg, my and m^). This is the usual MSSM set of input parameters, 
augmented by an additional CP-odd mass mA^^ and bn = the square of an angle parametrising Rp 

^ jl A V = implies there is no Rp violation at tree level in the -ino mass matrices. 
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violation. We define A2 to be the CP-odd scalar that becomes the z/ as (5r — * 0. Which of the Ai 
becomes the u is important, because we expect the Rp violating effects to go to zero as mj> 00 for 
all values of 6^. 

It can be seen from eqs. (|2^) and (^) that requiring the invariant S = v- [ISA ] ■ v/x^ tohe real, 
dictates that the mass splitting \m\_^ — m\J and 6r are not completely independent. Fixing the value 
of this mass splitting will give a maximum allowed value of from requiring that 5* be real; also the 
mass splitting must be greater than a certain value for a fixed value of 6r. 

In figure ^ we plot m^^ and as a function of 6r for values of niAi = 200, 300, 500, 1000 GeV 
with tan /? = 2, 10 and = 100 GeV. We observe the dependence of the maximum allowed value of 
6r on the mass splitting \m\_^ — m\J. As the mass splitting increases the maximum value of Sr also 
increases. Note that rrih^imh^) in figure 3 as — for tan/5 = 10(2) because this is the Higgs 

which becomes the sneutrino in this limit. We see that the lightest mass eigenvalue m/j^, decreases 
with 6fi for fixed CP-odd masses. On the other hand, m^^ increases as a function of 6r. The effect of 
6r on the heaviest eigenvalue m/13 is not very strong: we obtain m/13 ~ rriA^ for all the allowed values 
of 6r. 

Conversely, in fig. ^ we present the variation of m^^ with respect to for 6r = 0, 0.2, 0.5, 0.8, 
having fixed m^^ = 1 TeV, and for two values of tan (3 = 2, 10. For each value of 6^ there is a maximum 
allowed value of m^a- Recall that A2 is a sneutrino component when 6ji = 0. As m^a — > 0, so does 
ruhi because the lightest CP-even Higgs is the mode that becomes when 6r = 0. As increases, 
the mode "that would be the MSSM h ii 6r = 0" becomes the lightest Higgs and the plot flattens 
out. In the plot nih^ does not become exactly zero for 6r = 0, = due to one-loop corrections 
proportional to A'^ from the squark-quark sector. 

As mentioned above, the loop corrections induced by values A' ~ 1 are similar to those induced by 
the top Yukawa for the Higgs which couples to the up-type sector. The other i?p-violating coupling 
we have introduced in the calculation fii is always constrained by the neutrino mass to be sufficiently 
small that its contributions are negligible. 

Including i?-parity violation in the Higgs sector can be understood as having two effects on Higgs 
production and decay. It mixes the "Higgses" with the "sneutrino", and allows new decay modes for 
the Higgs/sneutrino decay products. There is of course no distinction between a Higgs and a sneutrino 
in the presence of i?-parity violating couplings; by "sneutrino" we here mean the CP-even and -odd 
mass eigenstates that become the sneutrino in the 6r ^ limit, and the "Higgs" is the CP-even mass 
eigenstate that would be the Higgs in the same limit. We define A2 to be the CP-odd scalar that 
becomes the uj as Sr 0. Mixing the Higgs with the sneutrino means that the eigenstates hi can all 
be produced via Z Zhi and Z hiAj, where i : 1..3 and j : 1..2. All of the hi can decay to bb, and 
to x^^j X^T ^iid XX if these decay modes are kinematically accessible. 

The neutralino can decay in the detector, via its production vertex (the neutralino becomes a 
neutrino and an off-shell Higgs, which can then decay to SM fermions). So unless 5r is uninterestingly 
small, the x^J^ and x^'x" should be visible. 

The new Rp violating decay modes hi — >• x'^t aiid hi x°z/ have been previously discussed ^ 



We plot the branching ratio to as a function of 5r in figure (|^). We assume in this plot that the 
decays to x+r and are kinematically forbidden. As expected from equation (p6D, the decay rate 
increases with 5r. The decrease at large 5r is a consequence of our parametrisation. 5r is sin^ of 
the angle between the vectors v and M^; as the angle increases to 7r/2 for fixed m^^ and m^g? \^u\ 
decreases. So the Rp violating mass term \Mu\^/5r. decreases. For larger values of tan/3 the decay 
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hi bb is dominant. 

Suppose now that the neutrahno is also heavier than hi, so only the Standard Model decays are 
available to the Higgs. The Rp violating couplings can still affect the production cross-section of the 
Higgs, and therefore the experimental lower limits on rrih- 

The production cross-section for Z — >■ Zh can be parametrised by = a{Z Zh)/a{Z Zh)sM- 
See, e.g., for experimental limits on In the MSSM, = sin^(/3 — a). It can be very small in 
our Rp violating model because it goes to zero a.s 6^ for the CP-even Higgs that becomes 9 in 
this limit. 

If is heavier than the CP-even Higgs which becomes h of the MSSM in the 6r ^ limit, then 
the CP-even Higgs corresponding to in the same limit will be heavier than (see figure In 
this case the Z Zhi vertex (equations ^ and ^9]) does not differ much from its MSSM value. We 
plot ^ as a function of 6r on the RHS of figure |^ for m^^ = 100 GeV. The present experimental lower 
limit on the Higgs mass for ^ ~ .8 is a few GeV below the ^ = 1 limit of 95.2 GeV ||56[| . 

Alternatively, if is light [], then can be very small. For instance on the LHS of figure ^ we 
plot ^ for the CP-even Higgs which becomes part of the sneutrino when 6^ 0. As expected, ^ is 
very small for small 6ji, because sneutrinos in the MSSM are pair-produced. 

Decreasing the ZZh vertex would decrease the experimental Higgs mass bound from this process; 



for ^ < .3, there is virtually no experimental lower limit ||56[. However, T{Z hA) increases as 
r{Z hZ) decreases, so there should still be a bound on m/j. The experimental lower limit on 
nih from Z hA is not trivial to determine, because the vertex and the two scalar masses are 
independent parameters. There are experimental limits in the MSSM |^^, but in this case the vertex 
and rrih determine m^. There are also bounds on sneutrino masses from Z uu* in models 
with trilinear Rp violation, but these assume that the CP-even and CP-odd sneutrino components are 
degenerate. The experimental lower limits on the Higgs masses in this model are therefore unclear, 
but likely to be lower than in the MSSM. 



7 Conclusions 

We have described the i?-parity violation induced by the additional soft mass terms in the scalar sector 
in terms of a basis-invariant quantity R (or Sr) . This eliminates the ambiguity usually present in these 
models when a specific Lagrangian basis for the hypercharge —1 doublets is chosen. We have analysed 
the effects of the i?p-violating couplings on the CP-even and CP-odd scalar masses to one-loop in 
a basis-invariant way. We have also calculated the Rp conserving and -Rp violating branching ratios 
of the lightest Higgs boson as a function of the basis-invariant quantity Sr. We have identified the 
regions of parameter space for which the decay modes of the Higgs boson are not those of the Standard 
Model Higgs. We have also calculated the production cross section as a function of Sr, and found that 
this can be strongly modified with respect to the Rp conserving case when the lightest Higgs boson 
is mostly "sneutrino-like" . The LEP lower bound on the Higgs mass in this model can therefore be 
lower than in the MSSM. 

''We assume that x^'^ decays are nonetheless kinematicaUy not allowed; the stau is therefore the LSP, but it decays 
so is not cosmologically a problem. 



17 



Acknowledgements 



S.D. and M.L. would like to thank J.R. Espinosa, P. Gambino, G. Ganis, G. Giudice and C. Wagner 
for useful discussions. N.R. wants to thank the CERN TH Division for kind hospitality. This work was 
supported in part by DGESIC under Grant No. PB97-1261, by DGICYT under contract PB95-1077 
and by EEC under the TMR contract ERBFMRX-CT96-0090. 



Appendix A: (z>) ^ basis 



We take the potential V to be the sum of equations (|T^) and ([15[) . We define v = \Jv1 + v\ = y/vjv^, 
and 

tan/3=-. (45) 

V 

We write the stop and sbottom masses as 



where for the stops 



Ml = ml + hlxV'^ + {9'-\9"){v'-x')l^ 



and 



For the sbottoms 



and 



= {Atx + htfi ■ v) 



Ml = ml + {X-vr/2-{g' + ^g''){v^-x')/^ 
Ml = ml + {\-vf/2-g'\v'-x')/l'2 



(46) 

(47) 
(48) 

(49) 

(50) 
(51) 

(52) 



The one loop minimization conditions can be expressed in terms of the CP odd mass matrix Mjj 
as in equations ( ^2] ) and (|23|): 

Muiv^ 
M A — 



X 








(53) 
(54) 



The CP-odd mass matrix is: 



Mud 

MuL 



( Mud MuL 
( Mdd Mrfi 



(55) 
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where the components are 



Muu = 

Mui = 



m: 



, a, mi COS 2/5 ^ 



+ 



b Ab^ 



We have defined 



m 



where Q is the renormahsation scale in the MS scheme, and 



A 

Db 



3 1 

3 1 
32^ 



with 



A 2 2 A 2 2 



(56) 
(57) 



(58) 
(59) 

(60) 
(61) 

(62) 



The CP-even scalar mass matrix is: 



M, 



uJ 



M 



7J 



2 ...2^ , /^,.2i_H< 



M„„ + m^sin^/5+<^— log ^ 



3/1,2 2A,Xa.rnl 



167r2 At ° m? 



log -^ + AtXt g{mt^,mij + {Xm ) X^g{mi^,mi^ 



t2 



vj 



Muj — rriz cos (3 sin (3 h 



167r2 



(^;,A^)(/.^A^)9^1og 



m 



b2 



167r2 



m. 



t2 



(63) 



M,^ + m|cos2/3^ + 



/^J + [^b{^K^i^)g{mf,^,mi^) 
3 



2 2 



ml 



A' 



A/A, 



XthtAtg{mi^ , m,-^ ) j [(64) 



+ 



167r^ Aft m~ 



167r2 

(AM^, + A,4)+ X,2<7(mft^,m^2) 



A'jA'} 



(65) 
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where 



g{mi,m2) 



IGtt^ {ml — 



mf — 7712 ^2 



(66) 



Appendix B: (z>) = 



This Appendix contains one-loop formulae for the minimisation conditions, the CP-odd mass matrix 
and the CP-cvcn mass matrix, in the basis where the sneutrino vev (u) = vl is zero at one loop. 
We define the up- type Higgs vev to be xl ^i ^-^^d the down- type vev to be v/\/2, so 

tan/3 = - (67) 

V 

and Xt = AtX + htfiv. 

In this basis, we can safely neglect the loop corrections due to hb and the soft trilinear coupling 
A'^ oc hh, since they are constrain to be small by the 6-quark mass rrib — —hbv/\/2 <^ v. If ml = 
[m|^](id,m| = and 771.3 = ^-re tree-level Higgs mass terms, /i = /iq {e = /ii) is the Rp conserving 
(violating) superpotential mass, and A' = A\^, then the minimisation conditions are 

ml — — mgtan/? — -m| cos2/3 -|- (68) 

ml = -777,3 cot /3 + ^mz cos 2f3 + 5ml (69) 

[m|]^^ = -BLt&n/3-hteDt^ - A'X'et&nPDh (70) 

where 

Sml = -htiiDt— , (71) 

V 

Sml = -^/i?[/(m,-J + /(7n,-J-2/(7n,)]-AA^-A'VA. (72) 

We have defined Dti and f{m) as in the previous appendix. 
The CP-odd scalar mass matrix elements are 

Muu = - cot (3{ml + htnAtDi) (73) 

Mud = ml + htnAtDt (74) 

MuL = BL + htAteDt + A'eWDb (75) 

Mm = -iQ.n(3{ml + htnAtDt) (76) 

= [mlU + hliieDt (77) 

2 

Mll = m\^+m\^ + {ml + ht^lAtDt)—— (78) 

Sill ■^A-' 
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with 

= - ^ { ("^Ai + "^A2 ) cos + 2hti^AtDt | 



2 



2 V ' ' COS2/3 
The CP-even scalar mass matrix M' is: 

= - cot p{ml + htiiAtDt) + ml sm^p + X'^e\^g{mi^,m-J 



+ 



167r2 



rriz 



2 



<<i = ^3 - ^ Sin 2/3 + /it/x A A + r^^r^X log + /itMt^t ^(^fi , ) 

t2 



/,3 \/3,,2, 



M'dd = - tan /3 {ml + htfiAfDt) + m| cos^ P + h^i/X^g {m^^ , m^^ ) 
M'dL = rnlL + hl^eDt + hlXl^eg{mi^,mi^) 

where and g{mi,m2) are as defined in the previous appendix. 

Appendix C: Some equations 

The normahsation factors for the basis-independent Higgs mixing angles, at tree level, are 

u 1 



det[N'] ^(dei N')2 + y'2' 

where N' = I - M', and 



For h — hi,h2 or h^, 



det N' = n^h~ ''^1(^1 cos^ /? + m\^ + m\J + sin^ f^'^\jn^A2 + '^'(^I — "i^l) 
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and 



y'2 = ml{SHan^P/il-5R)+2mlSsm^p + m%cos^psm^P) 

+m\jn\^ sin^ [3 cos^ [3 + 2m| cos^ (3 sin^ /^(m^^m^^ — m\m\){m\^ + m^^ — S/ cos^ /3) 
+m|cos2/?sin2/5[(mi^ +m^^ - Sf - 2m\^m%sm^ p - SHan'^(3/{l - 5r)] + 
m%S'^ cos2 /3sin2 - bn) - 2m| sin^ /^m^^^^^y^^ _ 

-2m2mi^mi^ sin^ + 2m| cos^ /?) - 2S'^m\ sin^ - 

(' 

For the CP-odd Higgses, the normahsation factor is 



[9 

[lo: 
[11 

[12 

[is: 



^9) 



n = ^ = (90) 

where N = I - M, = N -e ■ M„, a = either Ai or A2, and 

det N = - ml{Tn\^ + m^^ - 5*) + m\^m\^ sin^ /3 (91) 

and 

T/2 = m4^2^^^2^/^^ _ _ 2mlm\^ml^S sm"^ p + m\m% sin^/^cos^/? (92) 
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101.0 





Figure 3: and ruh^ as a function of for m^^ = 200, 300, 500, 1000 GeV and m^ij = 100 GeV. 
The input parameters for the loop contributions to the difference between the CP-even and CP-odd 
mass matrices are mg = 500 GeV, m^/ = tjid = 300 GeV, A = 200 GeV, and fij = (200,0). In the 
plots on the left, tan /? = 2; tan /3 = 10 for the plots on the right. 
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Figure 4: m/u as a function of ruA^ for 5r — 0, 0.2, 0.5, 0.8 and uiai — 1 TeV. The input parameters 
for the CP-even — CP-odd Higgs mass difference are rriQ = 500 GeV, mu = tud — 300 GeV, A — 200 
GeV, and /jLj = (200, 0). tan/3 = 2 on the left, and tan/? — 10 on the right. 



27 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 



Figure 5: The branching ratio for hi x°z/ as a function of 6r, for different values of and tan j3. We 
take ruAi = 500 GeV, = 120GeV , and the input parameters for the loops contributing to the CP- 
even — CP-odd Higgs mass difference are A = 0, fij = (200, 0), mq = 500 GeV, mu = mo = 300GeV. 
We take the total decay rate to be Ttot = ^{h bb) +T{h 
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Figure 6: The ratio ^ as a function of Sr, where ^ = gzzh/gz¥h ^he ratio of the ZZh vertex its 
value in the SM. When the hghtest CP-even scalar hi corresponds to the sneutrino in the Sr ^ 
limit, ^ can be small and goes to zero with Sr. If m^jl'^hich becomes ma when 6r = 0) is large, ^ is 
near its MSSM value. 
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